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SUMMARY 


Availability of a suitable propulsion system Is generally acknowledged to be a key 
requirement for the successful development of a new airplane. This paper describes the 
computational techniques utilized at Lewis Research Center to determine the optlmun pro- 
pulsion systems for future aircraft applications and to Identify system tradeoffs and 
technology requirements. 

Over the last five years, the NASA Lewis Research Center has obtained a greatly In- 
creased capability of performing detailed studies of engine cycles on the computer. Many 
more parameters can now be accounted for In the engine selection process. He can calcu- 
late cycle performance, engine weight, predict costs and account for Installation effects 
as opposed to fuel consumptlrn alone. Almost any conceivable turbine engine cycle can te 
studied since we do not rely on preconfigured simulation codes but can Input the engine 
cycle externally to the codes. Most of this capability has come through the Joint ef- 
forts of the Naval Air Development Center, The Boeing Company and NASA Lewis. 

These computer codes are: 

NNEP - a very general cycle analysis code that can assemble and arbitrary matrix of 
fans, turbines, ducts, shafts, etc.. Into a complete gas turbine engine and compute on- 
and off-deslgn thermodynamic performance 

HATE - a preliminary design procedure for calculating engine weight using the com- 
ponent characteristics determined by NNEP 

LIPCYC - a computer code presently being developed In conjunction with the Navy to 
calculate life cycle costa of engines based on the output from HATE 

INSTAL - a computer code presently being developed under contract to calculate in- 
stallation effects. Inlet performance and Inlet weight 

POD DRG - a table look-up program to calculate wave and friction drag of nacelles 

Examples will be giver, to Illustrate how these computer techniques can be applied to 
analyze and optimize propulsion system fuel consustptlon, weight and coat for representa- 
tive types of aircraft and missions. 

INTRODUCTION 

The airplanes, engines and missions of today are far more complicated than those of 
Just a few years ago- The ability to determine the optimum combination of airplane and 
engine Is of paramount Importance. But what Is the optimum combination. Is It the en- 
gine that bums the leas*, fuel?; Costs the least to operate?; Can mlrdmlze Installation 
penalties?; Minimizes fuel plus engine weight?; A combination of the above? 

Each alrplane/englne system probably has Its own criteria of optimization. It Is 
therefore necessary to develop the analytical tools capable of calculating all the fac- 
tors which enter Into the selection process. This paper discusses the computer tech- 
niques employed at the BABA Lewis Research Center to perform these calculations. The 
process by which almost any conceivable turbine engine can be evaluated as to fuel con- 
sumption, engine weight, coat and Installation effects Is described. Examples are shown 
as to the benefits of variable geometry and of the tradeoff of fuel burned versus engine 
weight. Future plans for further Improvements In the analytical modeling of engine sys- 
tems are also described. 


HANDKATCHINQ 

In order to deterailne engine operating characteristics at specified flight condi- 
tions, methods ware developed In the 19*0’s for superimposing engine component matching 
m«p« for simple engines such as turbojets and turboprops. These methods Involved labori- 
ous band calculations and performance map transformations to determine at what operating 
conditions of the engine components continuity of mass and energy, and mechanical speed 
relationships were satisfied. Needless to say, especially when methods were developed 
for a two spool engine, hours and even days were required to determine an operating line 
for an engine. A thorough discussion of these methods can be found In reference 1. Fig- 
ure 1 Illustrates the time frame and capabilities that existed. 
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EARLY COMPUTER BASED MATCHINQ CODES 

With th« advent of high speed computers, the task of matching of the engine compo- 
nents could not only be solved faster but more complex engines such as tuo-spool engines 
with a bypass flow (turbofan) could be simulated. Many companies, universities, and 
government Installations developed computerized methods. One of the earliest of these 
matching computer codes was called SMOTE and was developed at Wright Patterson Air Porce 
Base (ref. 2). SMOTE was capable of matching two spool turbofan engines. This capabil- 
ity was expanded by Che development of GEMENG and GENENQ II at NASA’s Lewis Research 
Center. The GENENG codes (refs. 3 and A) were capable of matching one, two or three 
spool engines with as many as three nozzles. Turbofans with booster or supercharger 
stages on the compressors could be simulated as well as aft-fan engines. GENENG served 
as the main simulation code In NASA and was adopted for use by over 20 Government agen- 
cies, companies, and universities. A version of GE'IEIIQ called DYNCEN was developed at 
Lewis to simulate transient behavior of turbofan engines (ref. 5) for use In control sys- 
tem studies. 


THE NAVY/NASA ENGINE PROGRAM (NNEP) 

Since 1973, the NASA Lewis Research Center has been conducting studies of advanced 
supersonic engines Including Variable Cycle Engines or VCE's. These engines take advan- 
tage of the use of variable geometry components and In-fllght flow switching capabilities 
such as from mixed flow to separate flow to attempt to deliver good engine performance at 
supersonic conditions as well as subsonic. By optimizing the exhaust profile during take- 
off significant decreases In Jet noise can also be achieved. It became apparent that 
GENENG and similar codes could not simulate some of the concepts coming out of the stud- 
ies. The new cycles did not fit Into any of the engine concepts already built Into the 
codes. 

Two options were available. A new specific code could be developed for each new 
engine concept, or a general code capable of simulating any engine could be developed. 

The second alternative was chosen as being more time efficient In the long run and more 
responsive to any Immediate need. We, therefore, decided to develop a new computer code 
In which an arbitrary engine configuration consisting of selected combinations of compo- 
nents could be described at Input time. It was also necessary to allow changes In engine 
configuration while running the code to simulate the operation of various VCE concepts. 
Furthermore, because of the large number of variables. It was highly desirable to opti- 
mize the settings of variable cotponents such as nozzle or turbine areas (e.g. , to mini- 
mize SFC for a given thrust). 

Contact with the Naval Air Development Center, Warminster, PA, revealed that they 
had a computer code, NEPCOMP (ref. 6), which already contalr.-d some of the features de- 
sired and whose structure was flexible enough to permit the additl'r. of others. This 
code lacked optimization capability and the ability to operate with "stacked" maps which 
would represent variable component performance. However, It already had the capability 
for processing arbitrary engine configurations. NASA-Lewls therefore contracted with the 
Naval Air Development Center for the Joint development of an Improved computer code. The 
objective of the Joint effort was to obtain a code capable of: simulating any turbine 

engine the user could conceive, simulating variable component performance, changing air- 
flow paths while running, and optimizing variable-geometry settings to minimize the spe- 
cific fuel consumption or maximize the thrust. 

An Interim version of this new code given the acronym NNEP (Navy NASA Engine Program) 
became operational In Hay of 197A and has been continuously refined since then to include 
all of the desired capabilities. 

NNEP contains almost all of the subroutines ^and Incorporates the philosophy of con- 
struction of NEPCOMP as described In reference 6. The major Improvements Incorporated In 
NNEP relative to NEPCOMP are the addition of: (Da performance optimization capability, 

(2) processing of stacked component maps for VCE operation, (3) nultlconflguratlon (modes) 
to simulate flowpath switching, (A) a computer generated engine configuration schematic, 
(5) throttle dependent inlet and boattail drag calculations, and (6) a simpler Input data 
format . 

As previously mentioned, the engine Is configured at Input time In running NNEP. 
First, the user draws a schematic of the engine he wishes to study, for example, a simple 
turbofan as shown In figure 2. Be assigns a flow station n’jmber 1 at the entrance to the 
Inlet and labels the inlet as component number 1. After this he is free to assign any 
ntinber at the other flow stations In the engine and to label each of the components with 
any component ntaber. One problem that does arise Is that Is not possible at all times 
to label the flow stations In accordance with the Aerospace Recommended Practice ARP 755A. 

The components that can be simulated In NNEP are as follows: 

Flow components - falling under this classification are 


(1) Inlets 

(2) Duets/burners 

(3) Compressors 
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(4) Turbines 

(5) Klxers 

( 6 ) Heat exchangers 

(7) Splitters 

(8) !Jazzles 

(9) Water Injectors 

Mechanical cotiponents - 

(1) Shafts 

( 2 ) Loads 

Control and ootlmlcatlon components - 

(1) Controls 

( 2 ) Optimization variables 

(3) Limit variables 

There Is a limit of a total of 60 components (Including all of the flow, aeohanlcal, 
control and optimization variables) allowed within the code. The maxlmua number of an; 
one type of flow or mechanical components Is 24 and the maximum number of controls * 
optimization variables Is 20. A K0H?IG Input card Is then generated by the user for each 
component as shown In figure 3- This figure Is for the compressor In figure 2. The com- 
ponent Is Identified as component number 4, that It Is a compressor and that Its primary 
upstream flow station number Is 4, there Is no secondary upstream flow; that the primary 
downstream flow station Is number 5; and the secondary dowr.stream flow station Is number 
13 (bleed flow). After all the components have been "configured," KHEP generates Its cwn 
flow path logic by Joining components by the station numbers. 

Each component has associated with It up to 15 required Inputs describing the com- 
ponent. These Inputs are usually design values such as pressure rise or map numbers 
corresponding to prestored performance maps for the component. An Illustration of the 
specifications for the compressor In figure 3 Is shown In figure 4. 

Control Information Is also entered as Input Identifying both the Independent and 
dependent variable as shown In figure 5. Optliolzatlon variables are entered similarly as 
shown In figure £. 

NMEP has proven to be a powerful analytical tool. Its primary purpose Is to gener- 
ate engine performance data for mission analysis studies. A typical use Is shown In fig- 
ure 7. This figure Illustrates the specific fuel consumption of a supersonic turbofan 
engine as a function of er.glne thrust when the supersonic engine Is operated at a subson- 
ic cruise condition Nach 0.9 at 11 000 meters (36 009 ft). Shown on the figure are three 
curves. The bottom curve represents the engine performance on an uninstalled basis, that 
Is, a pure thermodynamic cycle calculation. None of the variable-geometry features of 
»he engine have been utilized. 

An engine and Inlet .hlch »lz“^ for supcrs.,nl" cruize can suffer significant In- 
stallation losses subsonlcly. At the reduced power settings the inlet will oe capable 
of swallowing more air than the engine requires resulting In Inlet spillage drag. The 
boattall aft of the engine will not be filled with engine air resulting In additional 
drag. Installed performance for the fixed-geometry engine Is represented by the upper- 
most curve. As can be seen the difference between Installed and uninstalled perfomuice 
Increases as englr.e thrust Is reduced. The engine specific fuel consumption Increases 
rapidly at the lower power settings. 

The Introduction of variable geometry features Into the engine can greatly change 
the shape of the Installed performance curve. The performance of the engine with a vari- 
able geometry nozzle and variable area low pressure turbine Is shown on the remaining 
curve. The optimization capability of NNEP has been used to determine the optimum values 
of the two Independent variables. As can be seen, the curve la essentially flat. The 
components have varied to maintain as high an airflow as possible through the engine to 
reduce the spillage and boattall drag. NMEP has proven to be a very versatile engine 
cycle computer code and Is cow In use at approximately 30 government Installations, com- 
panies and universities. 


WEIGHT ANALTSIS OF TURBINE ENGINES - WATE 

With NMEP we are capable of simulating almost any turbine engine cycle the user can 
conceive cf. Being able to calculate engine perforunce and hence the fuel consumed on 
a mission Is an Important part of calculating the vehicle performance. It Is also neces- 
sary to be able to calculate the engine weight, length and diameter. The engine weight 
represents a significant part of the empty weight of an airplane. The length and dlBBe- 
ter of the engine are Important In calculating friction and boattall drags, la order to 
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develop the eapahlllty, NASA Levla awarded a contract to the Boeing Military Airplane 
Development Dlvlalon of the Boeing Company to develop an engine wel^t estimation code. 

The first version of this code WATE-1 (ref. 8) was completed In 1Q77. It used a 
preliminary design approach where stress level, maximum temperature, material, geometry, 
stage loading, hub-tip ratio and shaft mechanical overspeed are used to determine individ- 
ual component weights. The total engine weight was then calculated as the sura of the in- 
dividual components. The contract required that the code predict both Individual compo- 
nent and total engine weight within tlO percent accuracy. 

A relatively high level of detail was found necessary In order to obtain the re- 
quired accuracy. Component weight data for 39 different engines were used as a data base. 
This data base Is shown In the figure 8. The list of engines Includes military and com- 
mercial, turbofans and turbojets, augmented and dry, hardware engines and proposed en- 
gines, and supersonic and subsonic engines. 

WATE 1 was constructed to operate as an adjunct to NNEP. After running a cycle 
point on NNEP the thermodynamic properties were fed to the WATE-1 set of subroutines 
along with Inputs representing the design features of the components. The engine weight, 
length, and dlmenslcns were then calculated. At the same time, parts counts are gener- 
ated for the engine such as number of blades, sice of discs, etc. 

In 1978, NASA Lewis awarded a follow-on contract to the Boeing Company to extend the 
capabilities of WATE 1. This new version, WATE 2 (ref. 9) was completed In 1979 and has 
added many desirable features. Weight determination Is done for each component at Its 
critical operating point as follows: NNEP Is now used to "fly” the engine throughout the 

flight envelope of the aircraft and the maximum values of the flow, temperature, pressure 
and engine speed stored for use In sizing the components. Based upon these critical con- 
ditions, the weight Is determined. The capability to calculate the weight of radial flow 
components and of small engines was added In conjunction with a subcontract to the 
Garrett Dlvlalon of AlHesearch Manufacturing Company of Arizona. The engine center of 
gravity and moments of Inertia are also now calculated. 

ne accuracy of the code Is shown In figure 9. As can be seen, all of the engines 
fall within the ±10 percent band and. In most eases, approach ±5 percent or better espe- 
cially In terms of engine weight. 

WATE has bullt-ln default values for most of the Inputs. If the user does net enter 
values, these default values are automatically used. Many of these were used In the cal- 
culation of these weights. If more Information was available to us, especially In terms 
of geometry Inputs of the rotating components, these already small errors could probably 
be reduced even further. 

The cocblnatlcn of WATE and NNEP Is a very powerful analytical tool. As an example, 
a recent study considered the question of optimum cycle parameters for a duct burning 
turbofan for a supersonic cruise airplane (ref. 10). Some of the results of this study 
are duplicated here. The fuel mass and bare engine mass for 88 950 newton (30 000 lb) 
thrust engines flying 6440 kilometers (4000 mile) operating at Mach 3.4, 16460 meters 
(54 000 ft) Initial atltude are shown In figure 10. These masses are shown as functions 
of Bypass Ratio and Overall Pressure Ratio (OPR) with and without duct burning. The 
cycle analyst looking only at the fuel, mass In figure 10 would conclude that the optimum 
engine would operate dry and have an OPR of about 16 at a Bypass Ratio of 1.8 or more. 
However, when the mission analyst adds the fuel and engine masses as shown In figure 11, 
the optimum engine operates with the duetburner on, an OPR of 13 and a BPR of 0.8. 


LIFE CTCLS COSTING - COST/LIPCYC 

The question of cost Is entering more and more Into the selection process for opti- 
mum engines. The Initial cost Is not the only criteria for selection. Total life cycle 
cost Including maintenance, spares, operating costs, etc. must be considered for many 
applications. In order to develop the capability of calculating Life Cycle Cost, NASA 
Lewis contracted with the Naval Air Development Center (NADC) In 1978 to receive their 
costing model. NADC In turn subcontracted with Boeing to supply then with the production 
cost of the engine. As previously mentioned the weight code WATE calculates parts counts 
and weights as well as total component and total engine weight. These weights are trans- 
ferred to cost estimating routines which are based on correlations developed by NADC and 
Naval Air Systems Command. This procedure Is flow diagramed In figure 13. The correla- 
tion parameter Is based on a system of classifying materials by similarity of applica- 
tions In engines (ref. 11). In this procedure, materials used In Jet engines are placed 
In one of a total of six relative coat categories having to do with a combination of 
manufacturing cost and raw materials cost. Carbon steel and aluminum are assigned the 
lowest classification and used as a reference. High-strength hlgh-temperature nickel 
cobalt alloys which are costly and difficult to machine are placed In the highest (fifth) 
classification. Because of peculiar differences Ir. cost and machlnablllty, titanium 
alloys are assigned a separate (sixth) classification. Two Induces are developed for 
each material class, namely 

(1) Relative material cost 

(3) Relative machining cost 

The product of these two Induces Is called the "relative weighing factor." 
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In the cost estimation procedure, the estimated weight of each engine component la 
first converted to raw material weight. A raw material weight to finished material 
weight scaling factor, referred to as "Buy /Ply" ratio, has been estimated for each com- 
ponent for state-of-the-art and for advanced production methods. Haw material weight Is 
then multiplied by the relative weighing factor, and the sun of all such component prod- 
ucts Is formed. The summation for all engine components Is called the "Maurer factor" 

In honor of Its originator, R. J. Maurer. The production cost of the engine Is estimated 
by the linear correlation (fig. 13) between engine manufacturing cost and the Maurer fac- 
tor (ref. 11). 

This code Is Just becoming operational at .NASA Lewis and no results have as yet been 
generated except for Isolated check cases In which predicted costs have been compared to 
t.he actual and appear reasonable. A final report should be published sometime during the 
summer of 1979. 

Having determined the engine cost. It la now passible to determine the Life Cycle 
Cost based upon the .NADC Life Cyole Costing Model. The Interrelationship of COST with 
LIPCYC Is Shown In figure lU, NADC will supply the l.nputs and nodels to calculate all 
the parts of the pie other than manufacturing costs. It Is anticipated that this wcrk 
will be performed In the fall of 1979. 


IMPROVED i;:STALLATIO:i HFPECTS .MODEL - I.'ISTAL 

The previous example (fig. ■■) of varying engine flow by the use of variable geometry 
to reduce Installation effects showed the Importance of Inlet and nomole component per- 
formance, external as well as Internal. That figure was generated with a simplified 
model for Inlet and boattall drag that Is built Into the MMEP program. 

It was decided that NASA Lewis should obtain a more sophisticated nethod for these 
calculations. Consequently, a contract was awarded to Boeing Ir. IcmS to provide a bread 
subsidiary program for determining power-dependent Inlet and afterbody Installation ef- 
fects and also Inlet/nacelle weights (presently not In the *'ATE code). Nochle welg.hts 
are already calculated within WAVE. In addition to generating the component performance 
maps for NNEP, the code can be Interactive with the cycle and hence a tradeoff of Inlet, 
afterbody, and cycle can be utilized In the design process. 

The types of performance maps to be generated are shewn In figures 15 and 16 for In- 
lets and.noccles respectively. The necessary maps are obtained from either a data base 
or theoretical calculations (ref. 13). The data base contains perforttance data (usually 
experimental) for a spectrum of Inlet (axlsymetrlc, CD, pitot, mixed ocm-presslon) or "-cc- 
zle (axlsymetrlc , 2D, twin, etc.' types. A derivative procedure (ref- Ic) can be used to 
adjust the data base fco changes in design Maoh number, sldeplate shape, subsonic dif- 
fuser loss, cowl lip bluntness, takeoff door area, external cowl Initial angle, bleed 
system design, and bypass system exit design. Items not Included In the data are deter- 
mined analytically. .Nozcle/afterbody data are treated In a similar manner. The data 
base, being primarily experimental j offers Increased confidence In areas that are diffi- 
cult to treat theoretically such as viscous effects. 

After selecting the 'nlet sice or slcl.ng .Mach number, t.he Inlet and nozzle are 
matched to the .N.NEP cycle davo and the Installed performance calculated as well as the 
respective weights. Trade off studies can now be made of such effects as the best cen- 
blnatton of bypass and spillage for minimum specific fuel ccr.sumptlcn. 

The final report for the work being performed under this stuly contract Is scheduled 
to be published In the early fall of 1979- 


WAVE A!ID FRICTION D.RAG - PODDRG 

Under contract to NASA Langley Research Center, Rockwell International developed a 
method of evaluating the effects of nacelle shape on drag and weight cf a supersonic 
cruising aircraft (ref. H). Under this contract, Rockwell determined wave and friction 
drag Increments for a range of parametric shapes. As part of a follcw-on contract with 
NASA Lewis Research Center, Rockwell developed a computer code (PDDDRS) capable of Inter- 
rogating the data points generated under the previous contract In order to determine 
drag Increments for any nacelle shape of Interest (ref. 15). 

These nacelle Incremental drags are only applicable to the JLiSA arrow-wing superstn- 
le transport configuration (ref. 16). The program yields the Incremental wave and fric- 
tion drags of nacelles as functions of nacelle geometry variables and airplane Mach nijj- 
ber. The drag Increments are for the total vehicle relative to the vehicle with nacelles 
removed. That Is, all Interference effects with the airframe are accounted for. The 
nacelle shape parameters used as Inputs to the program are: 

(1) A^ Inlet capture area 

(2) Aj^j Nacelle maximum cross-sectional area 

(3) A^ Nozzle exit area (supersonic cruise position) 

(4) Distance from Inlet cowl leading edge to maximum crcas-sectlcnal area 



6 


(5) 1 Nacelle total lenftfc 

(6) Reference wlr^ area 

The output of this prcgrar. Includes, for the rs?tllt cf lr.'er*st: 

(1) The aforementioned Input data 

(2) Drag coefficients at Kach 1.2, rich and Input Mach number for fric- 

tion (CDFj, wave (CD’«3, and tctal ("i: irt.-a 

(3) The nondlmenslonal parameters of pcsltlcr. cf nixlnur. cross-sectional area 

(XaHAX/i)» noccle-to-capcure area rati -r nsslr.ur.-to-captur- area 

ration <Akax''“c 5* flr.er.esi ratlcr. 

In addition. Incremental drag coefficients of the reference airplane nacelle 
(ref. 16) are printed. 

Typical nacelle Incremental wave drag varlatlcrj are shown In figure 17. Mote that 
a properly shaped nacelle can produce a lower total airplane drag than that of the air- 
frame alone {LCZv * nacelle drac? calculated are fed Ir.te mission flight 

computer codes In evaluatl.-.g systems perforr.inoc-. 


CQKCL'JDINa REMARKS 

SASA lewis Research Center with a eomtlnstlor. of Ir-r.ouse, Joint, and contrasted ef- 
forts has beer, and Is continuing to develop the otran::* to determine the enelnee for 
opolmim mission performance. In the selettlcn process we can account for cycle perfor- 
mance, er.glne weight, life cycle costs, and installation effects. Future efforts will be 
directed towards Improvement In the current capatllltl-r , mainly in the areas of develop- 
ing better optimization methods to reduce eemputer time and analytl-ally determining per- 
formance maps for rotating machinery. For eaarple, we art about to award a contract for 
turbl.'.a Bsap generation, to be used for new cycles in which simple scaling of pre-eslitlng 
maps Is not sufficient. Ve believe that all of these r.c’.e.l' yl’l greatly reduce the ef- 
fort expended In performing clsslor. analysis by r.arrcwlr.r In mere oulekly on the er.glne 
cycles of greatest Interest. 
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1960 - • 2 SPOOL TURBOfAN (BYPASS ENGIND 

COMPUTERIZED MATCHING 

• TURBOJETS 

• 2 SPOOL TURBOFANS 

• 213 SPOOLS TURBOJETS, TURBOFANS, AFTFANS, ETC. 


1970 h 
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Figure L - Approximate history of methods of matching turbine 
engines. 






Figure 2 , • Simple 2 spool 





with cooled turbines. 















FLOW STATION NUMBERS 

Figure 3. - Define component type and location In (lowstream. 
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Figure 4 - Defining component ciaracteristics (for a conpressor). 
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Figures, - Owning controls. 
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Figure 6. - Defining optimization variables. 
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Figure 7. - Engine specific fuel consumption as a function 
of engine thrust. 
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Figure 8. - Data base engines. 
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Figure 11. - Sum of engine plus fuel mass. 88 950 N 
(20 000 lb) thrust engines. 
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Figure li - Overall program structure. 
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Figure a - Maurer factor correlation with cost 
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Figure 14. - Engine life cycle cost 
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Figure 15. - Format for inlet performance characteristics maps. 



% 


(s) DATA FORMAT FOR NOZZLE/AFTBOOV DRAG. 
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(bl DATA FORMAT FOR ROUND NOZZII GROSS THRUST 
COEFFICIENT. 



(d DATA FORMAT FOR NOZZLE GROSS THRUST CO- 
EFFICIENT FOR TWO-DIMENSIONAL NOZZLES. 


Rgure 16. - Format for nozzle/afUxxf/ drag and 
Cpg maps. 
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Figure 17. - Typical nacelle incremental wave 
drag variations with Mach number A. • 
2.79sqm(30sqft),//dc* 5.5. 


I 




% 



$ 


f 



End of Document 



